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Within massive gravity, we construct a gravity dual for insulator /metal phase transition and colos- 
sal magnetoresistance (CMR) effect found in some manganese oxides materials. In heavy graviton 
limit, a remarkable magnetic-field-sensitive DC resistivity peak appears at the Curie temperature, 
where an insulator/metal phase transition happens and the magnetoresistance is scaled with the 
square of field-induced magnetization. We find that metallic and insulating phases coexist below 
the Curie point and the relation with the electronic phase separation is discussed. 


I. INTRODUCTION 


In recent years, the holographic correspondence [1— 
4], relating a weak coupling gravitational theory in a 
(d + 1)-dimensional asymptotically anti-de Sitter (AdS) 
space-time to a d-dimensional strong coupling conformal 
field theory (CFT) in the AdS boundary, has been ex- 
tensively investigated and some remarkable progresses 
have been made in condensed physics systems [5-9]. 
For a recent review on the holographic superconduc- 
tor/superfluid models, please refer to Ref. [10]. Very re- 
cently, the present authors and their collaborators have 
realized the paramagnetism/ferromagnetism and param- 
agnetism/antiferromagnetism phase transitions in holo- 
graphic models by introducing a massive 2-form field in 
an AdS black brane background and some interesting 
magnetic properties of the models have been investigated 
in a series of papers [11-17]. In this Letter, we will pro- 
vide a new application of the holographic AdS/CFT cor- 
respondence by implementing the metal/insulator phase 
transition and the colossal magnetoresistance (CMR) ef- 
fect found in some manganese oxides materials in a holo- 
graphic model . 

Complex magnetic materials showing strong magne- 
toresistance have simultaneously been the focus of the 
attentions of the magnetic recording industry and the 
study of strongly correlated electron systems. Particu- 
larly, the study of the manganites such as Ay_,B,MnO3 
(A= La, Pr, Sm, etc. and B = Ca, Sr, Ba, Pb), which 
exhibit the “colossal” magnetoresistance effect, is among 
the main areas of research in strongly correlated elec- 
tron systems [18-24]. These materials show remark- 
able magnetoresistivity and an insulator (or semiconduc- 
tor)/metal phase transition associated with a paramag- 
netic/ferromagnetic phase transition, which has a com- 
pletely different physical origin from the “giant” mag- 
netoresistance observed in layered and clustered com- 
pounds. These materials are currently being inten- 
sively investigated by a sizable fraction of the condensed 
matter community, and its popularity is reaching the 
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level comparable to the high-temperature superconduct- 
ing cuprates. 

After great efforts in recent years, mainly through com- 
putational and mean field studies for realistic models, 
considerable progress has been achieved in understand- 
ing the curious properties of those compounds. However, 
a fully quantitative understanding of the CMR effect is 
still a challenge, much work remains to be carried out 
and it is the subject of current active research [25]. The 
holographic duality provides an alternative method for 
this type of strong correlated phenomena. In this paper, 
we will make a first attempt to build a holographic model 
to understand the CMR effect. 


II. HOLOGRAPHIC MODEL 


Before presenting our holographic model, let us make a 
brief analysis about how to build a holographic descrip- 
tion for such a phenomenon. Firstly, before ferromag- 
netic phase transition happens, CMR materials are in in- 
sulating phase where DC resistivity is finite and increases 
with decreasing the temperature. So the translation sym- 
metry is broken otherwise no scattering happens and DC 
resistivity is divergent. More important is that, more 
and more results from experiments show that the insu- 
lating phase in CMR is charge ordered, in which charges 
are localized and form inhomogeneous structures [21]. 
To realize the inhomogeneity for the CMR effect is a 
challenge both in condensed matter theory and holo- 
graphic description. Fortunately, just as pointed out re- 
cently in Ref. [26], momentum relaxation by breaking the 
translation invariance can be achieved by introducing a 
mass term of graviton in the bulk so that the macro- 
scopic DC resistance becomes finite. This provides us 
with a very simple holographic model to study macro- 
scopic DC resistance in some inhomogeneous materials 
without involving some complicated computations. Sec- 
ondly, in general, only breaking the translational invari- 
ance cannot lead to an insulating resistivity. Meffort and 
Horowitz [27] proposed a simple framework to have the 
insulating behavior in general relavity without breaking 
translational invariance, where a real scalar field is cou- 
pled to an U(1) gauge field. But Meffort-Horowitz model 
is only valid in the case of zero charge density. Thus 


a natural choice to build a holographic insulator model 
with finite charge density is to consider the Meffort- 
Horowitz model in a massive gravity theory. Finally, 
motivated by our previous work about DC resistivity in 
the paramagnetism/ferromagnetism phase transition in 
the probe limit [17], the model with massive 2-form field 
coupled with Maxwell field shows metallic ferromagnetic 
phase in low temperatures. This provides us with a mech- 
anism to describe metallic resistivity after the ferromag- 
netic phase transition. 

Based on these considerations, we present the model 
with the massive 2-form field-Maxwell-dilaton theory in 
a massive gravity with a negative cosmological constant. 
The action can be written as, 
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Here L is the AdS radius and G is the Newtonian gravi- 
tational constant. Without loss of generality, we can set 
L = 1/16rG = 1. à,J,a,8 and go are all model pa- 
rameters. R is scalar curvature, A, is the U(1) gauge 
field with field strength Fy, = 2V,,A,). wy is a dila- 
ton field with squared mass m3. My, is a 2-form field 
with squared mass m? and a nonlinear potential V(M). 
dM is the exterior differential of 2-form field M,, and 
(dM)? = 9V),M,7jV"M"".Lmg is the mass term of 
graviton, where matrix K is defined in terms of the dy- 
namical metric g,, and a reference background metric 


fuv by 
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This is a special case of the formulation of the massive 
gravity theory [28, 29] presented in [30]. As shown in 
Refs. [26, 31], there is a position-dependent mass of the 
gravitational perturbations, 

m3(r) = —28—a/r. (3) 
Ref. [26] shows that, in the leading level of lattice, the 
graviton mass term is proportional to the square of mod- 
ulating amplitude multiplied by wave vector in the dual 
boundary lattice. This gives a physical meaning of gravi- 
ton mass in holographic model, i.e., ms describes the 
strength of inhomogeneity in the dual boundary theory. 
This consequence plays a key role in our holographic 
model. In this paper, we will pay attention to two lim- 
its, i.e., weak inhomogeneous and strong inhomogeneous 
cases, which correspond to the cases of m2 < 1 and 
m2 > 1, respectively. These two cases admit us to ob- 
tain DC resistivity in a simple manner. We will see later, 


only the case with large value of ms can give rise to the 
typical CMR effect, which agrees with the fact that there 
is a strong inhomogeneity with CMR effect in the man- 
ganite. 

The parameter À can be understood as the coupling 
strength between the polarization field M,,, and back- 
ground Maxwell field strength. In the effective action of 
string tehory, the value of the dilaton coupling parame- 
ter is taken to be gj = 1. The case with go = v3 corre- 
sponds to the 4-dimensional action by dimensionally re- 
ducing from the 5-dimensional Kaluza-Klein theory [32]. 
Nonetheless, it is helpful to set go as an arbitrary positive 
constant here so that we can see what the role of dila- 
ton field plays in the model. In addition, the choice of 
nonlinear potential V (M) is not unique. What we need 
is that there is a critical temperature, below which the 
xy-component of M,,, can condense. In this paper, we 
take the potential as follows, 


V(M) = (“Mu M")? = (M A M)P. (4) 


Here * is the Hodge-star operator. We choose this form 
just for simplicity. To be stable for both the bulk and 
boundary theories when Y% = 0, we require ma > 0 for 
all r [33, 34]. Following [34], we take the degenerate ref- 
erence background with fr = fyy = 1 and other compo- 
nents are vanishing. Although the reference background 
is degenerate, Ref. [34] showed that this two-parameters 
massive theory is ghost free for the case with 8 mass 
term, but it has not yet been proven for the case with a 
mass term. Therefore in this paper, we will take a = 0 
in order to avoid possible problems with causality. Note 
that this reference background leads to that there is a 
conserved energy but no conserved momentum current 
in the boundary theory. 


II. BACKGROUND EQUATIONS AND 
PERTURBATIONS 


Now we are in the position to calculate the DC resistiv- 
ity from our holographic model (1). For this, we assume 
the dynamic metric has following form, 


dr? 
ref 
where f and y are two functions of r. Suppose the so- 
lution has a horizon at rp, the associated temperature 


then is T = r? f'e7x/2/4r. We take the following self- 
consistent ansatz for matter fields, 


Ay, = o(r)dt + Brdy, % = V(r) 
Mav = —p(r)dt A dr + p(r)dax ^ dy. 
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Here B is a constant magnetic field and it will be viewed 
as the external magnetic field in the boundary field the- 
ory. Put the ansatz and the dynamic matric in (5) into 
the action (1), we can get a set of ordinary differential 


equations. Near the boundary r — oo, the equations give 
the following asymptotic solutions for matter fields, 
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where 6; = y1+4m? and 62 = /9+ 4m, u is the 
chemical potentail, ø is the charge density and p+ and 
w+ are all constants. We impose the regular condi- 
tions at the horizon and Dirichlet and source free con- 
ditions for matter fields at the boundary of r > o, i.e., 
=f, pyro? = A and p4 = 0. Without loss of gen- 
erality, we can set u = 1. Note that nontrivial solution 
for ~ always exists when go and F, are both nonzero. 

Following Ref. [17], we need that J < 0 so that p 
can spontaneously condense below a critical temperature 
when B = 0. We also need that 6, > 1 and ô> < 3, 
otherwise, the nonlinear terms of p and w will play more 
and more important roles when r > oo, which will break 
the asymptotic AdS4 geometry of space-time and lead to 
the instability of the dual theory in the UV regiem. 

Now let us study how DC conductivity is influenced by 
temperature and external magnetic field in this model. 
Because of the planar symmetry in the boundary, the 
conductivity is isotropic. By the AdS/CFT dictionray, 
we need only turn on a perturbation of gauge field such 
as 0A, along the x-direction with ingoing boundary con- 
dition at the horizon. 


In the case of weak inhomogeneity ms <1, the main 
part of DC resistivity is very simple. In that case, the 
graviton mass is very small, gravitation fluctuations suf- 
fer from smaller scattering than others. So the DC con- 
ductivity is dominated by the background geometry. In 
other words, we can neglect the fluctuations of matter 
fields when we compute the DC resistivity. Then follow- 
ing Refs. [26, 31], we can find that DC resistivity can be 
expressed as 


Riigne = mg /[W?(1 — A? /4m3)] + O(m). (8) 


The more interesting case is the strong inhomogeneous 
limit, i.e., me > 1, which is the case we are really in- 
terested in this paper. In such a limit, graviton has 
very heavy mass so that it is in fact very hard to be 
excited by fluctuations. Heavy mass term suppress the 
fluctuation of gravity so that we can fix the background 
geometry. In this case, the main part of DC resistiv- 
ity can be obtained by just considering the fluctuations 
of matter fields. In the DC limit of w — 0, we need 
only consider three perturbations 6A, = €a,(r)e~‘“* and 
M;; = €C;;(r)e~“' with ingoing conditions at the hori- 
zon, where (i, j) = (r, x), (t,y). At the linear order of e, 


we have following equations in the low frequency limit, 
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where a prime stands for the derivative with respect to 
r, and w,p and p determined by the background equa- 
tions, O(w) is the terms with order of w, and the other 
equations of gravity parts and matter parts are order of 
O(w) , all of which can be neglected when w —> 0. 

At the boundary r — oo with p,p — 0, we have the 
following asymptotic solution for az, 


Ay 


Ay = Az4+ + par (10) 
Following the AdS/CFT dictionary, we can obtain that 
electric current is (J) = as- and the DC resistivity is 
given by R = lim,,-,9 iwa,4/(J). 

In fact, with the “membrane paradigm” of black hole, 
we can directly obtain the DC conductivity from Eqs. (9) 
using the method proposed by Iqbal and Liu in Ref. [35]. 
To see this, we first note that, 


lim r? f(a, — \°C,2/4) = —(1 — A? 4m? I. (1) 


Too 


Eq. (9c) shows that this quantity is conserved along the 
direction r when w — 0. At the horizon, using Eqs. (9) 
and the fact that Cy, is regular at the horizon when w + 
0, we have, 


(7) = iwa, (Tp) e- 2G0%bo dm? 
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(12) 
Here Xo, Yo, Po and po are the initial values of Y, p and p 
at the horizon, respectively. In the low frequency limit, 
Eq. (9c) implies that the electric field is constant, i.e., 
lim,=,r, Gz (1) = az+. So we obtain the DC resistivity in 
heavy graviton limit as, 
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As a self-consistent check, we can take A = 0. In that 
case there are only dilaton and Maxwell field. Then we 
have Rieavy = e790 + O(1/m2), which agrees with the 
exact result R7' = Z(Wo) +H? /m? given in Ref. [31] with 
daliton coupling Z(po) = e790. From the expression 
for DC resistivity, we see that when Rheavy ~ m2, the 
heavy graviton limit is broken. In that case, the fluctua- 
tions of metric have to be taken into account. 


IV. METAL/INSULATOR PHASE TRANSITION 
AND MAGNETORESISTANCE IN STRONG 
INHOMOGENEITY 


The physical phase in different temperature depends 
on the model parameters. As a typical case , we fix 
m? = Vane = —2, m? = 40,g9 = 1 and à = 3/4 to 
compute the DC resistivity at different temperature and 
small external magnetic field numerically. We first solve 
the equations of motion from the action and then calcu- 
late the DC resistivity numerically. All the results are 
shown in Fig. 1. 

In Fig. 1(a), we plot the DC resistivity at zero mag- 
netic field with different dilaton source A. There is a 
critical A, ~ 0.103 (we scan A from 5 to -1, numeri- 
cal precision restricts our ability to check wider region). 
When A > A,, there is an insulator/metal phase tran- 
sition. The resistivity shows an insulator’s behavior de- 
scribed by dilaton field when T > Te. When the temper- 
ature is lowered to the Curie temperature Tc, p begins 
to condense spontaneously and a ferromagnetic phase 
transition happens. Below and near To, the resistivity 
decreases when temperature is lowered, which shows a 
metal’s behavior. Though the behavior of DC resistivity 
is transformed into metallic from insulating, the insulat- 
ing phase described by dilaton field coexists with ferro- 
magnetic metallic phase in the sample. Numerical results 
show that two different electronic phases can coexist be- 
low the Curie temperature. There is a distinct peak at 
the temperature where spontaneous magnetization be- 
gins to appear and an insulator/metal phase transition 
happens there. This is just one of characteristic prop- 
erties of CMR materials in manganese oxides. When 
A < Aea, DC resistivity shows a metallic behavior in the 
whole temperature (up to the region that numerical com- 
putations can be done), though there is still a saltation 
at To. What’s more, when a small magnetic field B is 
turned on, we find that the resistivity is very sensitive to 
the external magnetic field (see Fig. 1(b)). 

Here we emphasize that the heavy graviton limit plays 
a very important role. Just as mentioned above, for 
the light graviton case, the behavior of DC resistivity 
is dominated by fluctuation of graviton, then no such 
metal/insulator phase transition or CMR effects appear. 
This is agreement with the fact that the CMR effect is 
due to the fact that two electronic phases mix with each 
other and form inhomogeneity in nm scale [19, 21, 36]. 

It is interesting to compare our holographic results 
with some experimental data of CMR materials. The 
resistivity of a typical CMR material La,_,Sr,MnO3 
is shown in Fig. 1(c) and (d). We see that our holo- 
graphic model gives qualitatively similar results when 
x > 0.1, which is a powerful evidence to support that 
this holographic model is a suitable one to describe the 
CMR. effect. Furthermore, in La,_,Sr;MnO3, when 
x > Ze ~& 0.2, the peak of DC resistivity disappears, 
which is very similar to the case when daliton source 
A < Ae. So this similarity gives us an explanation for A 
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FIG. 1. (a): The behavior of the DC resistivity vs temper- 


ature in different A with B = 0 and J = —2. (b): The 
DC resistivity vs temperature in different external magnetic 
field. Here Bo/T ~ 5.5 x 1074, and J = —1/3. (c) and 
(d): The DC resistivity for Lai_2zSrzMnO3 as a function of 
temperature for different doping x and magnetic field. The 
experimental data are from Ref. [18]. 


at the dual boundary field, i.e., dilaton field describes the 
effects of doping. In addition, we can find from Eq. (13) 


the following scaling relation for magnetoresistance (MR) 
in the case of weak magnetic field and T > T ; 


MR(B) = 1 — R(B)/R(B = 0) œ pĝ x B?. (14) 


Note that in the region of T > Tce, the system is in the 
paramagnetic phase. Therefore the magnetic moment 
N is proportional to B in the weak field case. Then 
Eq. (14) tells us that MR(B) « N?. This result is in 
complete agreement with the experimental data of CMR 
materials [18]. 


V. DISCUSSION 


In this paper we have present a gravity dual to the 
metal/insulator phase transition and found that the 
model can describe the CMR effect in some manganese 
oxides materials. The behavior of DC resistivity is in 
complete agreement with experimental data. The model 
provides a new example to apply the AdS/CFT corre- 
spondence to condensed matter systems. As the first at- 
tempt to describe the CMR effect in a holographic setup, 
more aspects of the model should be further studied. The 
first one is to study the behavior of DC resistivity in the 
case with arbitrary graviton mass. For this, we need 
to consider the perturbations of the gravitational back- 
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